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ABSTRACT 


The  authors  have  developed  an  interactive  graphics  system  for 
the  correction  and  reduction  of  airborne  laser  profiles  of  sea  ice. 

The  system  is  implemented  on  a Vector  General  renewable  display  system 
connected  to  a PPP11/45  computer,  but  minor  modifications  would  make 
it  suitable  for  any  computer  with  interactive  graphics.  The  laser 
profile  is  displayed  in  sections  and  the  human  operator  is  armed  with 
routines  for  the  removal  of  noise  spikes,  phase  jumps  and  aircraft 
altitude  variations.  Completely  automatic  profile  correction  techniques 
have  been  found  to  b<  inadequate  or  unworkable,  especially  in  heavily 
ridged  ice.  Although  designed  specifically  to  treat  the  output  of  the 
Geodolite  laser  prof i lometer,  the  system  can  be  applied  to  many  other 
data  correction  problems  in  geophysics  and  other  sciences  where  long 
time  series  of  data  are  obtained  with  periodic  noise  spikes  or  a long- 
period  drift  or  modulation  which  requires  removal. 


. - — rj 

W.ite  Section  tr 
N i,  «i  Section  O 

I vne  r. 

1 n a n 


I r • in 


INTRODUCTION 


■ 


' 

t 


A knowledge  of  the  surface  topography  of  sea  lee  Is  ol  groat  import- 
ance in  the  modelling  of  po tar  ocean  dynamics.  The  drag  coefficient 
between  wind  and  ice  Is  determined  by  the  surface  roughness  (Bankc  et 
a 1 , 1U7(,)  , and  since  most  ol  the  heat  flow  from  ocean  to  atmosphere 
occurs  through  thin  ice  (Maykut,  ll'7f>)  it  is  valuable  to  measure  the 
distribution  of  surface  elevations  as  a guide  to  the  distribution  of 
ice  thicknesses.  The  airborne  laser  prof i lomet er  is  the  most  effective 
instrument  toi  swiftly  generating  such  data  over  long  stretches  of  the 
ice  cover.  It  produces  a profile  of  surface  elevation  along  a Hack 
directly  below  t tie  aircraft,  using  a beam  whose  footprint  is  essentially 
a moving  point  and  which  is  therefore  capable  of  very  high  resolution. 

The  earliest  and  most  extensive  laser  surveys  ol  the  Arctic  Ocean 
wore  carried  out  by  the  Ub  Naval  Oceanographic  Office  (Ketchum,  1971; 

Welsh  and  Tucker,  1971;  Tooma  and  Tucker,  1 7 3 ) and  furnished  data 
which  were  used  to  test  theories  of  pressure  ridge  height  and  spacing 
distributions  (Hiblet  et_  al_,  197.’,  1974;  Diachok,  1975)  and  to  delineate 
regions  ol.  the  Arctic  Ocean  with  differing  degrees  of  roughness  (Tucker 
and  Westhall,  1973;  Hibler,  1975).  The  same  agency  lias  carried  out 
laser  profiling  in  the  Beaufort  Sea  as  a contribut  ion  to  t lie  A 1 IV)  EX 
project,  and  co-operative  profiling  with  a submarine  in  the  East  Greenland 
current  (K07.0  and  Diachok,  1973).  Laser  profiling  in  the  Arctic  Ocean 
has  also  been  carried  out  by  an  Argus  aircraft,  of  Department  of  National 
Defence,  Ottawa  (Dunbar  and  Lowry,  1974;  Wadhams,  197b).  As  well  as 
measuring  ridging,  it  is  also  possible  to  use  a laser  profi lomet or  to 
measure  attenuation  rates  ot  ocean  waves  in  the  marginal  ice  none 
(Wadhams,  1975).  To  date  no  laser  profiling  has  been  done  in  the 
Antarctic,  although,  sea  ice  characteristics  there  dilfei  significantly 
from  those  in  the  Arctic  and  a profiling  experiment  would  be  ot  great 
interest.  The  most  recent  lasei  profiling  experiment  in  the  Arctic 
Ocean  has  been  a eo-operat  ive  project  in  October  197b  involving  an 
Argus  aircraft  ami  a British  nuclear  submarine  following  identical 
tracks  at  approximately  the  same  time  (Wadhams,  1‘’77;  Wadhams  and 
Lowry,  1977),  thus  enabling  the  statistical  relationship  between  surface 
and  bottom  topography  to  be  derived.  The  interactive  graphics  system 
described  in  the  present  report  was  developed  to  analyse  the  data  from 
this  experiment. 

bo  tar  all  airborne  laser  surveys  have  employed  the  same  instrument, 
the  bpect t a-i’hysics  Geodoltte  *A.  This  uses  a continuous  lle-Ne  laser 
beam  with  a wavelength  of  b32.D  nr.’  and  an  output  power  of  ap\  roximately 
.’5  mW.  The  outgoing  beam  Is  sinusoidally  modulated  by  an  electro- 

optical  modulator,  and  after  reflection  from  the  target  the  beam  is  [ j 

received  by  a Schmidt -Cassegrain  telescope  ot  o..'  m diameter  which  is 
coaxial  with  the  laser.  The  phase  shift  between  the  modulations  on 
the  transmitted  and  received  beams  is  measured,  giving  the  range  to  the 
target.  A set  ol  modulation  frequencies  is  available,  usually  differing 
by  factors  ot  It’.  This  enables  the  Inset  to  be  used  ns  an  accurate 
1 1 timet  or:  the  lowest  modulation  f regucncy  gives  a full  scale  deflection 
(O  to  .’it  in  phase)  tor  an  altitude  change  ot  It’  leet,  and  the  higher 
t 1 equestci  give  1 1 '** , h \ It’'  anti  10  t eet  lull  scale  deflection. 


By  switching  frequencies  stepwise,  each  significant  figure  of  the  altitude 
in  feet  can  be  obtained.  When  used  to  obtain  sea  ice  terrain  profiles 
only  the  higher  frequencies  are  used,  giving  better  resolution. 

Unfortunately,  the  aircraft  normally  porpoises  with  an  amplitude  exceeding 
lo  feet,  and  ice  ridges  can  also  exceed  10  feet  in  elevation,  so  that 
a laser  profile  on  the  10  foot  scale  has  frequent  "phase  jumps".  A 
phase  jump  occurs  when  the  range  crosses  one  end  of  the  10  ft  interval, 
ie  the  phase  crosses  0°  or  360°.  When  this  happens  the  electronics 
automatically  put  in  a 360°  jump,  positive  or  negative,  in  order  to  put 
the  profile  back  on  scale. 

Figure  1 (top)  shows  a typical  raw  profile  produced  by  this  type 
of  instrument.  There  are  four  components  in  the  profile:- 

1.  The  actual  sea  ice  elevation. 

2.  A profile  of  aircraft  porpoising,  in  which  both  altitude  and 
attitude  variations  combine  to  give  a smoothly-varying  range  oscillat- 
ion of  wavelength  about  2-4  km.  (ie  Figure  1,  centre). 

3.  Phase  jumps,  which  are  theoretically  instantaneous  and  of  identical 
amplitude.  In  practice,  because  of  the  recording  system  used 
(magnetic  tape)  and  imperfections  in  the  prof ilometer , the  jumps 
have  a finite  rise  time  of  up  to  0.05  s and  are  of  varying  amplitude. 

The  profilometer  used  for  the  1976  experiment  was  modified  to  give 
a 20  ft  full  scale  deflection  on  its  most  sensitive  scale,  reducing 
the  frequency  of  phase  jumps,  and  to  give  a jump  of  180  instead  of 
360°.  This  restores  the  profile  to  the  centre  of  the  scale  instead 
of  to  the  opposite  extremity,  reducing  the  likelihood  of  the  profile 
immediately  jumping  back  by  the  same  360  . 

4.  Noise  spikes  and  "drop-outs"  (none  visible  on  fig.  1).  A drop- 
out is  a brief  positive-and-negative  spike  caused  by  a momentary  loss 
of  phase  lock.  The  loss  may  be  due  to  the  aircraft  crossing  an  ice- 
water  boundary,  or  else  experiencing  an  exceptionally  sudden  range 
change  (eg  an  overhanging  ice  block  in  a ridge) . Other  spikes  may 
be  caused  by  interference  from  other  equipment  on  the  aircraft.  Longer 
periods  of  noise  can  occur  due  to  phase  lock  failure  over  thin 
cloud  or  undercast  mist. 

The  purpose  of  a data  processing  system  is  first  to  remove  (3) 
and  (4),  leaving  a profile  as  in  fig.  1 (centre),  and  then  to  remove 
aircraft  porpoising  as  accurately  as  possible  so  as  to  leave  a true 
profile  of  ice  surface  elevation  as  in  fig.  1 (bottom) . Computer 
programs  can  be  written  to  automatically  remove  phase  jumps,  but  the 
process  is  very  unsatisfactory  since  the  profiles  must  always  be 
examined  afterwards  and  sections  with  imperfect  jump  removal  or  with 
noise  spikes  present  must  be  discarded  or  further  treated.  Similarly, 
there  is  no  fully  satisfactory  way  of  automatically  removing  aircraft 
motion.  A simple  filter  will  not  work,  since  the  data  do  not  represent 
a gaussian  random  process,  the  deviations  from  the  porpoising  being  all 
in  the  positive  direction.  A three-stage  filter  was  described  by 
Hibler  (1972)  to  overcome  the  last  problem.  The  method  used  is  to  connect 
up  the  local  minima,  thus  underlining  the  data.  This  underlining  is 
then  filtered  and  subtracted  from  the  corrected  data  to  yield  an  ice 
profile.  The  main  disadvantages  of  this  procedure  are  that  it  does  not 
allow  for  the  existence  of  lec.ds,  and  again  that  it  does  not  allow  for  interactive 
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viewing.  Thus  it  must  be  coupled  with  an  automatic  phase  jump/noise 
spike  removal  program,  which  works  only  in  regions  of  good  quality  data 
(ie,  normally  lightly  ridged  areas) . 

The  system  described  in  this  report  is  based  on  that  developed  at 
Defence  Research  Establishment  Ottawa  by  Lowry  and  Brochu  (1976)  and 
uses  a human  operator  to  do  the  pattern  recognition  for  both  errors  and 
aircraft  porpoising  removal,  leaving  the  computer  to  manipulate  the 
data  and  do  all  the  necessary  computations.  Data  is  presented  on  a 
television-type  display  and  the  operator  indicates  the  position  and 
nature  of  errors  with  a light  pen.  When  the  errors  are  removed  to  the 
operator's  satisfaction  the  data  is  underlined  by  the  operator,  again 
using  a light  pen.  The  underlined  data  is  then  filtered  and  subtracted 
from  the  corrected  data,  leaving  a profile  of  pure  ice  elevation,  free  of 
aircraft  motion.  The  corrected  profile  is  written  to  magnetic  tape  so 
as  to  be  available  for  statistical  analysis  on  the  same  computer  or  on 
a separate  general-purpose  computer. 


2.  DESCRIPTION  OF  THE  PROGRAM 


The  system  has  been  implemented  on  a PDP11/45  computer,  operating 
under  DOS,  with  a disc  drive,  two  magnetic  tape  drives  and  a Vector 
General  visual  display  unit  with  light  pen  and  keyboard.  The  program 
was  written  mainly  in  Fortran. 

The  program  sequentially  displays  blocks  of  data  from  an  input 
tape;  these  are  interactively  corrected  by  the  operator  and  the 
result  displayed  for  inspection.  If  satisfactory,  the  blocks  are 
output  onto  a second  magnetic  tape.  During  this  process,  sufficient 
data  is  kept  on  disc  so  that  the  operator  can  stop  at  any  time  and  pick 
up  later  where  he  left  off,  and  so  that  there  is  minimum  loss  should 
there  be  a system  crash. 

The  main  program  keeps  track  of  the  block  numbers,  controls  the 
disc  and  tape  files,  and  calls  the  correction  subroutines.  After 
asking  the  operator  whether  the  data  is  new  or  old  (ie  system  being 
restarted  after  a break)  the  program  reads  in  the  relevant  block  of 
data  into  the  array  ICE,  which  also  holds  subsequent  corrected  versions 
of  the  data.  The  program  then  hands  this  over  to  the  subroutine 
CORREK  for  the  correction  of  phase  jumps  and  noise  spikes,  and  then  to 
the  subroutine  UNDERL  for  the  drawing  of  the  aircraft  motion.  This 
second  routine  produces  an  array  PTS  holding  the  drawn  underline. 
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This  underline,  which  consists  of  line  segments,  is  filtered  by 
a Hamming  filter,  which  essentially  produces  a weighted  running  mean 
of  200  points.  However,  to  avoid  edge  effects  of  the  filter  at  the 
beginning  and  end  of  a block,  the  array  PTS  must  be  preceded  by  at 
least  lOO  points  of  the  underline  from  the  previous  block  (stored  in 
array  BSTR } , and  must  be  followed  by  the  first  lOO  points  of  the  subsequent 
block  (stored  in  array  ESTR) . The  very  first  and  last  blocks  on  the 
input  tape  must  therefore  have  end  ramps  constructed  for  them.  For 
this  reason  the  filtering  (subroutine  HAMFIL)  must  take  place  a block 
behind  the  correction  and  underlining.  Thus  for  new  data  the  sequence 
of  events  is:-  I 

— process  blocks  1 and  2 

— filter  and  display  block  1 

— process  olock  3 

— filter  and  display  block  2 etc. 

When  a block  is  filtered  and  displayed,  the  operator  may  reject  it, 
in  which  case  this  block  is  treated  as  a first  block  of  data  and  the 
program  restarted,  or  he  may  accept  it,  in  which  case  the  smoothed 
underline  is  subtracted  and  the  result  similarly  displayed  for  approval. 
When  finally  approved,  this  block  is  written  onto  output  tape. 


It  should  be  noted  that  the  program  has  been  linked  with  three 
standard  libraries  of  routines  held  permanently  on  the  Cambridge 
PDP11/45  computer:- 


1. 


2. 


3. 


A standard  Fortran  subroutine  library 
Magnetic  tape  routines 


MTREAD 

MTWRT 

mtwtm 

MTSKIP 

MTRWND 


reads  a block 
writes  a block 
writes  a tape  mark 
skips  blocks 
rewinds 


Graphics  routines 
NITDEV 
NITBUF 
BCNPIC 
ENDPIC 


initialise  device 
initialise  a picture  buffer 

for  creating  a picture  segment 


DELPIC 

WKiTOL 

INWAIT 

REATOL 

LINI 

LINE 

LINIR 

WINDW 


for  deleting  a picture  segment 
for  positioning  tracking  square 

for  inputting  on  visual  display  unit/Vector  General 
keyboard 

for  finding  position  of  tracking  square 

absolute  and  relative  line 
segment  drawing 

for  setting  up  a display  window 


The  subroutine  CORREK  has  two  main  functions.  Firstly  it  seeks 
large  jumps  between  consecutive  points.  It  indicates  these  in  turn, 
asking  the  operator  to  say 
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i)  that  it  is  a phase  jump  or  noise  spike,  in  which  case  it  is  removed,  or 
ii)  that  the  jump  is  genuine  data,  in  which  case  it  is  kept. 

After  this  the  routine  waits  for  commands  from  the  operator  to  remove 
remaining  errors,  redraw,  reject  and  start  again,  or  finally  accept  the 
modif ications  (see  OPERATOR'S  GUIDE). 

The  subroutine  UNDER!.  waits  for  commands  to  insert  or  delete  a 
line  segment,  alter  the  starting  point  of  the  line,  or  finally  accept 
the  underline. 

A basic  block  diagram  for  the  main  program  is  shown  in  Pig  2. 

Block  diagrams  for  the  subroutines  will  be  found  in  section  4. 

Some  modif ications  to  the  simplest  form  of  the  program  were  necessary 
because  ol  space  and  speed  limitations  on  the  PDP11/45.  The  program  fills 
most  of  the  available  45K  of  store,  so  space  saving  procedures  were 
employed  which  makes  the  program  (section  4)  slightly  difficult  to 
follow.  Further  KAMFIL  as  written  in  FORTRAN  takes  90  secs  for  a filter, 
and  so  was  replaced  by  an  equivalent  routine  written  in  Assembler 
code,  which  requires  only  IS  secs. 

The  Hamming  Kilter  was  chosen  tor  its  low  side  lobes  in  the 
Fourier  transform  domain.  It  should  be  pointed  out  however,  that  a 
simple  top  hat  filter  could  be  used.  Its  first  side  lobe  is  only 
13dB  down,  but  in  many  cases  this  is  entirely  adequate.  (If  the  data 
are  to  be  subject  to  Fourier  analysis,  however,  the  side  lobes  become 
more  important) . A top  hat  filter,  once  filled,  will  filter  any 
length  of  data  block,  with  only  a few  simple  calculations,  thus  taking 
much  less  time. 


AN  OPERATOR'S  GUIDE  TO  THE  BY  STEM 

This  section  refers  specifically  to  the  operating  system  of  the 
PDP11/45  - Vector  General.  It  is  included  to  illustrate  the  relative 
ease  of  carrying  out  those  operations  on  an  interactive  graphics  system. 
The  program  is  set  up  on  the  user's  disc;  the  data  is  in  digital  form 
on  magnetic  tape;  the  output  is  written  to  a second  magnetic  tape. 
Normally  Geodolite  data  is  recorded  on  analogue  magnetic  tape,  so  that 
A-to-D  conversion  is  necessary.  The  I97(i  data  was  recorded  thus  on  a 
Tandberg  Berios  hX1  KM  instrumentation  tape  recot d -i  with  flutter  compen- 
sation, and  digitally  converted  using  an  AiX'l-D  analogue-to-digital 
converter  built  into  the  PDP11  '4S  computer.  Digital  data  points  were 
O.Ol  s apart,  corresponding  to  about  0.9  m of  track;  an  analogue 
low-pass  filter  was  included  in  the  A-to-D  converter  to  remove  all 
energy  above  the  Ny^vtsf  frequency. 


INITIALISE 


Figure  2. 


../NEW 

"\DATA? 


READ  ANCII 
FROV  DISC  I 


PROCESS 
NEXT  BLOCK 


KILTER  PREVIOUS 
BLOCK 


SUBTRACT 

UNDERLINE 


WRITE  BLOCK 
\ TO  TAPE  / 


CLOSE  DOWN 


Simplified  block  diagram  of  the  main  program 
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Initialisation 


The  operator  loads  the  input  tape  (data  tape)  on  unit  O and  the 
output  tape  on  unit  1.  Me  loads  the  disc  on  the  disc  drive  and  initialises 
the  PDP11  with  the  sequence 

HALT 

LOAD  ADDRESS  173100 
LOAD  ADDRESS  177406 
START 

The  Vector  General  display  is  switched  on  and  adjusted  for  high  gain 
(since  the  light  pen  responds  more  readily  at  high  intensity) . 

When  the  PDP11  is  initialised  the  operator  logs  in  on  the  Teletype 
keyboard  using  his  password  (eg  LO  61,166)  and  assigns  the  Vector 
General  to  output  stream  6 

AS  VG : , 8 } 

then  starts  the  program  running  with 
RU  Gj 

where  G is  the  name  of  the  program. 

Operations 

1.  Upon  initialisation  of  the  proqram  the  operator  is  presented  with 
the  question 

NEW  DATA  OK  OLD?  (N  OR  0) 

in  the  top  left  hand  corner  of  the  screen  (where  all  comments  or  questions 
from  the  program  to  the  operator  appear) . lie  answers  on  the  Vector 
General  keyboard  with 

01  if  he  wishes  to  continue  with  a dataset  begun  in  an 
earlier  session,  or 

Nj  if  he  wishes  to  start  a new  dataset. 

For  new  data  the  program  sets  up  its  own  ramp  function  to  fill  up  one 
end  of  the  Hamming  filter.  Otherwise  it  searches  for  the  last  block  to 
be  treated  and  uses  the  last  part  of  that  block  to  fill  up  the  filter. 

2.  The  program  now  presents  a block  of  2.000  joints  (about  2 km  of  data) 
on  the  screen,  scaled  sc  as  to  fill  the  active  area,  with  a y=0  line 
included  (although  this  may  be  in  an  arbitrary  position  due  to  wrap- 
around) and  with  a vertical  bar  in  the  bottom  left  corner  representing 
the  scale  of  200  digitising  units  (about  0.66  m) . If  the  profile  has 
any  line  segments  for  which 

yin)  - y(n-l)  > 200 

the  trucking  square  will  position  itself  at  the  bottom  of  the  screen 
directly  below  the  jump  in  question,  starting  with  the  jump  nearest  the 
start  of  the  block.  The  operator  now  has  the  option  of  accepting  the 
program's  judgement,  ie  of  having  the  jump  removed  on  the  grounds  that  it 
is  a phase  jump  or  a noise  spike.  To  do  this  he  just  types 
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It  ho  thinks,  however,  that  the  jump  represents  real  data  (part  of  a 
ridge)  ot  that  it  can  be  better  removed  by  the  "0"  facility  (see  - below) 
ho  cat:  reject  the  program's  decision  by  typing 


N } 

The  tracking  square  now  moves  on  to  the  next  jump  until  all  jumps  are 
exhausted.  An  instantaneous  phase  jump  lias  an  amplitude  of  10  feet 
t~”lJ  digitising  units);  we  use  .'00  as  the  criterion  because  jumps 
seldom  occur  instantaneously. 

3.  If  there  are  now  critical  jumps  in  the  block,  or  when  the  operator 
has  dealt  with  all  those  that  exist,  the  progiam  presents  on  tin*  screen 

TRY  G R X OR  0 nnn 

where  nnn  is  the  number  of  the  block  (for  information  purposes).  Ii  the 
operator  types 

R ,) 

he  cat;  see  the  results  of  his  operations  in  stage  (j)  since  this  redraws 
the  block  in  its  updated  state,  le  with  the  results  ot  all  manipulations 
to  date.  R can  be  used  at  any  point  during  this  stage  of  processing  to 
examine  the  latest  state  ot  the  data. 

The  "G"  facility  deals  with  any  undesirable  features  of  tin*  data  that 
still  remain.  It  replaces  a short  length  of  data  with  a horizontal 
straight  line  segment.  It  can  therefore  be  used  to  draw  a line  to 
replace  a noise  spike,  or  else  to  take  out.  a phase  jump  which  occupies 
a finite  number  of  data  points  instead  of  being  instantaneous.  l'o  use 
this  facility  the  operator  uses  the  light,  pen  to  move  the  tracking 
square  to  the  lei t side  ot  the  proposed  line  segment  and  types 

G } 

He  then  moves  to  the  light  side  ot  the  proposed  segment  and  types 

For  typical  applications,  see  EXAMPLES. 

It  the  operatoi  is  not  satisfied  with  the  u It  ei at  ions  that  he  has 
made  he  can  type 

x: 

which  redraws  the  block  in  its  original  unmodified  state,  rejecting  all 
oi  tlie  operator’s  manipulations  to  date. 

When  the  operatoi  is  finally  satisfied  that  the  block  not*  consists 
only  of  real  data,  with  spikes  and  jumps  removed,  and  subject  only  to 
aircraft  mot  ion  (which  should  appear  as  a smoothly  trending  curve)  tie 
types 

Phe  program  new  proceeds  to  the  next  stage. 


4.  It  redraws  the  block  once  more  and  presents  the  message 


UNDERLINE 

It  also  presents  the  tracking  square  at  the  extreme  left  side  of  the 
screen  on  the  y-value  that  it  had  reached  while  underlining  the  preceding 
block.  If  the  operator  is  dissatisfied  with  this  as  a starting  position 
for  his  underlining,  he  moves  the  square  to  a better  position  at  the 
start  of  the  block  and  types 

B} 

Otherwise  he  simply  moves  the  square  along  (using  the  light  pen) , stopping 
at  points  which  he  feels  are  representative  of  the  "level  ice  surface" 
and  typing 

) 

In  this  way  he  constructs  a piecewise  straight  line  profile  of  the  surface 
of  level  ice.  To  recognise  this  surface  requires  judgement,  experience 
and  sometimes  imagination;  a prior  study  of  some  vertical  aerial  photographs 
of  the  Arctic  ice  cover  is  an  advantage.  The  main  problem  is  not  to  be 
misled  by  very  smooth  ice  (with  a "roughness"  which  is  only  the  noise 
level  of  the  original  tape  recorder) . This  is  often  thin  polynya  ice, 
and  the  operator  should  try  to  keep  a short  distance  above  such  ice, 
following  the  level  of  smooth  floe  ice  (which  usually  has  greater 
roughness) . Narrow  cracks  can  also  be  mistaken  for  level  ice  segments. 

The  aim  is  to  achieve  consistency;  ie  to  follow  either  the  level  of 
polynya  ice  throughout  (for  open  icefields)  or_  to  follow  the  level  of 
smooth  thick  ice  throughout,  which  is  more  appropriate  for  the  central 
Arctic  where  polynyas  may  be  far  apart. 

If  the  operator  is  dissatisfied  with  his  latest  line  segment  he 
can  type 

D i) 

which  deletes  the  last  segment  drawn.  Repeated  use  will  delete  part  or 
J whole  of  his  line.  If  by  mistake  he  does  not  proceed  monotonically  in 

the  positive  x-direction  he  will  receive  the  message 

MOVE  ALONG 

and  his  last  segment  will  be  disregarded.  If  he  exceeds  the  maximum 
allowable  number  of  segments  in  a line  (25)  the  program  sends 

DELETE  SOME 

and  the  operator  must  delete  some  of  his  points  and  redraw  the  rest  of 
the  line  with  points  further  apart. 

After  taking  his  line  to  the  extreme  right  of  the  block  (with  the 
help  of  a chart  recorder  printout,  if  available,  to  determine  the 
continuing  trend) , the  operator  types 

0.) 

which  sets  the  filter  routine  in  action. 


5.  The  program  now  filters  the  preceding  block's  underline,  with  the 
aid  of  the  first  few  points  in  the  current  block  to  fill  up  the  Hamming 
i liter.  When  the  filtering  is  complete  (about  15  seconds)  the  program 
presents  the  preceding  block  together  with  its  piecewise  line  filtered 
into  a smooth  curve,  and  asks 

OK?  Y OK  N 

The  operator  examines  the  smoothed  underlining  to  see  whether  it  appears 
to  be  a true  representation  of  the  level  ice  surface.  If  ho  thinks  that 
it  is,  he  types 

Y 3 

otherwise  N J 

A rejection  takes  the  program  back  to  stage  (4)  for  the  preceding  block; 
it  presents  this  block  again  for  underlining.  The  formerly  current 
block  is  lost  and  must  bo  processed  again  after  the  preceding  block. 

On  acceptance,  the  program  subtracts  the  smoothed  piecewise  profile 
from  the  laser  data  to  obtain  a final  laser  profile,  allegedly  a true 
representation  of  the  ice  surface.  It  presents  this  again  to  the 
operator  with  the  message 

OK?  (Y  OK  N) 

Again  the  operator  has  the  option  of  rejecting  it,  with  the  same  consequences 
as  before.  If  he  accepts  the  final  version,  it  is  written  onto  the  output 
tape  and  the  program  asks 

A REST?  (Y  OR  N) 

If  the  operator  wishes  to  wind  up  the  session  he  types 
Y) 

and  the  program  rewinds  both  tapes,  having  stored  the  current  operating 
data.  The  operator  now  types 

fO 

on  the  Teletype  keyboard,  which  stops  the  program  so  that  it  is  safe  for 
him  to  remove  tapes  and  disc. 

If  he  wishes  to  process  another  block,  he  types 

N}  on  the  Vector  General  keyboard 

whereupon  the  next  block  of  raw  data  is  read  in  from  the  input  tape  and 
presented  to  him. 


It 


4.  II.UISTRATED  EXAMPLES  OF  THE  USE  OF  THE  SYSTEM 


! 


I 


I 


In  this  section  we  show  several  photographs  of  the  Vector  General 
display  screen  which  demonstrate  aspects  of  the  operation.  Figure  3 
shows  the  calibration  signal  used  to  eliminate  errors  from  the  tape 
recorder  and  digitisation.  The  profilometer  output  could  be  held  at 
any  one  of  three  DC  values  corresponding  to  +180  offset,  zero,  and 
-180°.  At  least  5 seconds  of  each  of  the  three  signals  were  recorded 
at  the  beginning  of  each  tape  and  whenever  adjustments  were  made  to  the 
sensors.  The  left  hand  side  of  the  trace  shown  in  fig  3 has  a 0°  t.o  -180° 
calibration  signal.  The  average  value  of  each  offset  was  calculated 
from  a print-out  of  the  screen  display,  and  Used  to  establish  the 
relationship  between  the  unit  ot  digitisation  in  the  computer  and  the 
actual  elevation  perceived  by  the  profilometer  (the  absolute  accuracy 
of  the  profilometer  is  dependent  on  the  stability  of  d quartz  oscillator 
kept  in  a temperature  controlled  oven,  and  is  of  the  order  of  l part 
in  10s).  Figure  3 also  gives  an  impression  of  the  noise  level  introduced 
into  the  data  by  the  recorder  and  digitiser;  this  is  just  the  thickness  of 
the 'line"  representing  the  PC  signal. 

Figure  4 shows  a typical  section  of  heavily  ridged  ice  from  the  1976 
co-operative  profiling  experiment.  The  vertical  bar  at  the  lower  left 
corner  is  included  b^  the  pr  gram  to  give  the  operator  an  idea  ot  the 
vertical  scale,  since  the  computer  always  re-formats  the  profile  to 
occupy  the  whole  height  of  the  screen.  The  bar  is  200  units  (0.86  m) 
high.  The  tracking  square  has  placed  . t sol t automatically  beneath  a 
spot  where  there  is  a jump  of  more  than  200  units,  in  this  case  it  has 
clearly  identified  t he  position  ol  u l.'o1  phase  shift. 

In  figui  .>  the  opeiator  has  accepted  the  program’s  judgement  and 
by  typing  3 has  removed  the  jump.  The  jump  was  a simple  one,  and  its 
removal  reveals  the  presence  of  a high  pressure  ridge. 

There  is  a second  phase  jump  on  the  far  right  of  the  section.  This 
has  a finite  rise  time  and  was  not  found  by  the  program.  The  operator 
moves  the  tracking  square  to  the  beginning  and  end  of  the  jump,  and  by 
using  the  "G"  c mm and  reduces  the  tump  to  a short  horizontal  line  segment. 
Figure  6 shows  the  result  ot  this  operation. 

Since  the  raw  data  have  no  furtlui  phas.  umps  or  spikes  the  operator 
now  underlines  it  with  : rt  line  ;egment s joining  apparent  minima  or 
stretches  of  smooth  ice  (figure  7).  Figure  S shows  the  result  of  filtering 
this  piecewise  line  into  a s" ooth  euivv.  The  operator  accepts  the  result 
and  subtracts  the  smoothed  1 i ; t t >.  t prof  ile.  Figure  9 shows  the 
final  resu 1 t . 
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Figure  3 


Figure  4 


Figure  5 
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Figure  6 
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Figure  7 


Figure  0 
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MAIN  PFOGRAM  FCF  THE  CORRECTION  OF  LASER  PROFILES  OF  SEA  ICE 
COMMON  /A/W( A) /B/TC< 2) 

INTEGER  I CE( 2000) / PTSC 2000) / PI CBUF( 4200) / BSTR(  1 0 1 ) / ESTR(  101) 
INTEGER  ANCIH  105)/ V/ OFSET 

EC'l' I VALENCE  ( N/  ANCIL(  1 ) ) / ( OFSET/  ANCIL(  2)  ) / < LPT/ ANC I L ( 3)  ) 

EGl'I  VALENCE(M/  ANCILt  4)  ) / ( BSTR(  1 ) / ANCIL(  5)  ) / < BSTP(  1 0 1 ) / PTS ( 1 ) ) 
EQUI VALENCEC FT  SC  2000) / ESTRC  1 ) ) 


INITIALIZATION 


DEFINE  FILE  1 ( 1 / 1 05/ U/ NAN) 
DEFINE  FILE  2( 2/ 2000/ U/ NPTS ) 
DEFINE.  FILE  3(  2/  2000/  0/  N I CE> 
CALL  NI TDEV( 3 ) 

CALL  NITBUFCPICBUF/ 4200) 

CALL  MTRWND(O) 

CALL  MTRWND(l) 

CALL  HAMTAB 


C********************************** 

C STARTING  AFRESH  OR  CONTINUING  ? 


C*************1 


M = 0 

1 CALL  BGNPI C(  3) 

WHI TEC  8/  1 00) 

100  FORMAT ( 1H \,  ' NEW  DATA  OR  OLT?  CN  OR  0)  ’ ) 
CALL  ENT PI C 
READC8/101)  V 


101  FORMAT ( A 1 ) 

CALL  DEL PI C( 3) 

IF<  V. EC. ' N* ) GOTO  2 


RFADC 1*1)  ANCIL 

CALL  MTSKIPC 1/M/ IEFR/ IEOF) 

GOTO  3 

C STARTING  AFRESH 


2 CALL  BGNPI C( 3) 

VRI TEC  8/  106) 

106  FORMATC  1H  1/ ' STARTING  AT  E-LOCK  ') 
CALL  ENDPIC 
READC  8/  105)  N 
105  F0PMATCI4) 

CALL  DELFICC3) 

OFSE  T = 0 


LFT=0 

2 1 CALL  TAPINC I CE/N/ I EOF) 

CALL  CORREKC ICE/OFSFT/N) 
CALL  PNDE.RLC  PTS/  LTT) 

CALL  HAM1 GNC BSTR) 

K=MOLC  N/ 2) ♦ 1 
VRI TEC l’l)  ANCIL 
VRITECD’K)  PTS 
WRI TEC  3*  K)  ICE 


C»*******»************************************************************ 

C PROCESS  NEXT  BLOCK 

U •******»•***«**.«•»**•***»»***»****•****•*•  ******♦»•«  «*•«****••  *.«.*** 


n»s  page  is 
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3 N*N*  1 

CALL  TAriNCI  CE,N,  I EOF) 

I FC  I LOF. EC. 0)  GOTO  A 
CALL  HAMENDCESTR) 

GOTO  5 

4 DO  4 1 K*  I, 2000 

41  I CEC  K) * I CE(K) ♦OFSET 

CALL  COHREKC ICE, OFSET, N> 
CALL  DNDERLC PTS, LTT) 
K=MODCN, 2) ♦ 1 
WRITECP'K)  FTS 
WRI  TEC  3*  K)  ICE 
CALL  HAM COFC  TTS, ESTR) 

5 K = N-  1 
K=MODCK,  2) ♦ 1 
REA DC  2*  K > PTS 
READC 3* K>  ICE 
CALL  HAMFILCDSTP) 

CALL  VINDOVC I CE, 1/ 2000) 

CALL  DI  SPLAC  ICE,  1,2000,  2,  1) 
CALL  Dl SPLAC PTS/ 1/ 2000/ 2/ 2) 
C******************************** 
C IS  THE  DNDEHLINt  OK  ? 


to  CALL  BGNPICC3) 

WRI TF  C 8/  107 ) 

107  FORMATC  1H 1/ * OK?  Y OK  N *) 
CALL  ENDPIC 
READC  8/  101)  V 
CALL  DEL PI CC1) 

CALL  DELPI CC  2) 

CALL  DELTI CC  3) 

CALL  PELTICCOQ) 

IFCV.EO. ' N * ) GOTO  3 
DO  to  1 K= 1/2000 
to  1 I CEC  K>  = I CEC  K > - TTS C K ) 

CALL.  WINDOWC  I CE/  1/ 2000) 

CALL  DI  Sn.AC  ICE/  1/  2000/  1/  1 ) 


IS  THE  RESDLT  GOOD  E.NOl’GH  TO  GO  ON  TAPE  ? 


CALL  DGNPICC3) 

WRI TECH/ 107) 

CALL  ENDPIC 

7 READC8/101)  V 

CALL  PELTI CC  1 ) 

CALL  CELPI CC  3) 

CALL  PELTICC^*)) 

IFCV.EO. ' N * ) GOTO  8 

CALL  MTVRTC 1/ I CE/ 2000/ I ERR) 

M = M ♦ 1 

WRI TEC  1 * 1 ) ANCIL 
I FC I EOF.NE.O)  GOTO  11 
CALL  BGNTI CC  3) 

WRI  TEC  8/  103) 

103  FORMATC 1H 1/ ' A REST?  CY  OH  N) • ) 
CALL  F.NDriC 

10  RE  ADC  8 / 101)  V 

IFCV.NF.’Y*  • AN  D.V.NE. *N* ) GOTO  10 
CALL  DELPICC3) 

IFCV.EO. 'N* ) GOTO  3 
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C SUBROUTINE  F OH  REATING  A BLOCK  FROM  THE  INPUT  TATE 


SUBROUTINE  TAPIN(  I CF.,  N,  I EOF) 

INTEGER  ICE(l) 

TATA  NE/I/ 

IF(N.EC.NB)  GOTO  1 
K=N-NB 

CALL  MTSKIPCO,K, 1 ERR/  I EOF) 

I F ( I ERR.FG.O)  GOTO  1 I 
WRI TE<  6,  1 00)  N 

FORMAT ( * ERROR  ON  SKIPPING  TO  BLOCK', 15) 
STOP 

I F < I EOF • NF • 0 ) RETURN 
NB  = N + I 

CALL  MTREAPC  0#  ICE,  2000,  I ERR,  I EOF) 

DO  2 1=1,2000 
I CE  ( I ) = - I CF  < I ) 

IFC I FRR.FC.O)  RETURN 
WRITEtb, 101)  N 

FORMAT ( ' ERROR  ON  READING  FLOCK', 15) 

STOP 

END 


SUBROUTINE  T0  CORRECT  FATA  FOR  TM ASF  SHIFTS  ANP  NOISE  STIKES 


SUBROUTINE  CORHEK( I CE, OFSET,N) 

C OMM ON / A/W ( 4 ) /P/TC(  2) 

PI  MENS  I ON  X( 2>,Y( 2) 

INTEGER  1 ED<  2), ICE ( 1 ), OFSET,L,H 
PAT A 1 ED/ 2,-1/ 

KA=OFSET 

CALL  WINPOW< ICE, 1,2000) 

CALL  DI SPLA( ICE, 1, 2000, 1, 1) 

DO  44  1 = 1,  1399 
J»ICE< !♦ 1 > -I CE ( I ) 

IF(O.LT. 200.ANP. J.GT.-200)  GOTO  44 
TC<  1 ) = ( I -W<  1 ) ) /< W< 2) - W ( 1 ) ) 

TC<  2) =W( 3) 

CALL  WRI  TOL<  20  1 , 0,  0,  TC,  2) 

K=INWAIT(  - 1 .,  I ED,  I PA,  1) 

I F ( I PA • NF . ' ' ) GOTO  44 

JJ=I ♦ 1 


THIS  PAGE  IS  BF.ST  QUALITY  rRACriCABUj 
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C J CM 


DC  23  K = J J * 2000 
23  ICECK)=ICECK)-J 

CF5ET=0FSET-J 
44  CONTI  NUE 

CALL  EGNPI CC 3) 

OKI  TC  8*  101  ) N 

101  FORMAT!  1H  1*  * TRY  G H X OR  OS  10X^14) 

CALL  ENDFIC 

I I*INWAIT(- 1 . * I ED# IDA*  1) 

IFC I CA.NE. *G* ) GOTO  2 
CALL  KEATOL<  201*0*0*  TC*  2) 

XC  1 ) = TCC  1 )*< W<  2) -W(  1 ) ) *WC  1 ) 

Y<  I)=TCC2)*CW(4)-WC3))*WC3) 

K 1 = I F I X ( X < 1 ) ) 

I = INWAI T< - 1 .* I ED) 

CALL  REATOL< 201* 0* 0* TC* 2) 

X C 2) = TCC  1 )*C  VC  2) -W(  1 ) ) ♦ WC  1) 

YC 2) = TCC 2)*(W( 4) -WC 3) )+WC 3) 

K 2-  IF IX(X(  2)  ) 

I V=I FIXCYC  1 ) -YC  2) ) 

CO  11  KK  = K 1 * K2 

II  I CEC KK) =1 CE(K1 ) 

K2=K2+ 1 

TO  12  K K = K 2 * 2000 
12  I CECKK)=ICECKK)+IV 

OFF ET-OFFET+ I V 
GOTO  1 

IFC  IDA. ME.  *H*  .AND. IDA. NE.  'X'  .AND. IDA. NE.  '0’  ) GOTO  1 
IFC I DA. ME. *X'  ) GOTO  22 
CALL  DELPICC3) 

CALL  TAFINC I CE*N* I EOF) 

OF  SET=KA 
DO  21  KK= 1*2000 

1 I C EC  KK ) = I CEC  KK ) ♦ OF  SET 

2 CALL  DELPICC99) 

CALL  WINDOW! ICE*  1 * 2000 ) 

CALL  DELPI CC 1 ) 

CALL  DI SPLAC I CE*  1 * 2000*  1*1) 

IFC  IDA. EG).  'X*  )GOTO  45 
IFC I PA. EQ. * O' ) RETURN 
GOTO  1 
END 


C SUBROUTINE  TO  UNDERLINE  THE  DATA 
C ****'*****♦****** 

SUBROUTINE  UNDERLC PTS* LPT) 
COMMON  /A/W( 4) /P/TCC 2) 

INTEGER  I EDC 2) * PTSC  1 ) * I TMPC  1 00) 
TATA  I ED/2* - 1 /MAX/48/ 

CALL  DELPI CC 3) 

CALL  EGNPI CC  3) 

WRI TEC8* 100) 

100  FORMATC  1H  I *’ UNDERLINE* ) 

CALL  ENDPI C 
CALL  EGNPI CC  2) 

CALL  ENDPI C 

6 1 TMPC  1 ) = 1 FIXCVJC  1 ) ) 

I TMPC  2) =LPT 
K = 3 

CALL  DRAW C I TMP*  K ) 
i I a I N VAI  TC  “ 1 • * I ED*  I PA*  1 ) 


■■*  + + * t * * * * * 
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I Ft  I r/>  - r U - * CJ  * 3 GCi'0  4 
1 Ft  I DA. ME. * • ) GOTO  2 

IF(K.LT.MAX)  GOTO  II 
CALL  DELPI Ct 3) 

CALL  I3GNPI  C(  3) 

VRITEtS, 101) 

101  FOKriATC  1H \,  ' DELETE  SOME') 

CALL  ENDPIC 

GOTO  13 

11  CALL  REATOLt  20 1 1 0< 0#  TO*  2) 

ITMPtK) = I FIX( TCt  I )*(VK  2) -Wt  1 > )f Wt  1 ) ) 

I TMPt Kf 1 ) =1 F1X( TC( 2) * < W< A) -W(  3)  ) f Wt  3)  ) 

I  Ft  I TMPtK)  .GT.  I TMPCK-2)  f I > GOTO  12 
CALL  DE.LPI C( 3) 

CALL  DCNrI C<  3) 

WRITECa, 102) 

102  FORM AT ( 1 M 1 * 'MOVE  ALONG') 

CALL  ENCPIC 

GOTO  1 

1 2 K-K+  2 

CALL  PRAWCITMP/K) 

13  GOTO  1 

2 I Ft  I DA.NE. ’ D* ) GOTO  3 

IFCK.E0.3)  GOTO  1 

!<  = K-  2 

CALL  DRAW < I TMP*  K) 

21  GOTO  1 

3 I Ft  1 DA.NE. ' f ) GOTO  35 
CALL  REATOLt  201*  Oj  0»  TC.*  2) 

LFT  = I FIXtTCt2)*tWt4)-W(  3 ) > ♦ Wt  3)  ) 

GOTO  6 

35  CALL  rELPlCt3) 

CALL  EGNPI Ct  3) 

W H 1 TEtS, 103) 

103  F ORMAT  t 1H  1 < ' TRY  [-<  I > C<  OP  O') 

CALL  ENDPIC 

GOTO  1 

k CALL  DELPI Ct 1 ) 

CALL  DELPI Ct 2) 

CALL  DELPI Ct 3) 

CALL  EF.LPI  Of  99  ) 

TCt 1 ) =0. 

TCt  2) =0. 

CALL  WHI TCLt  20 1 » 0* 0/ TC>  2) 

1 TM  POO  =2000 
I TM Pt K+  1 > = I TMPt  K- 1 ) 

K = K-  2 

DO  4 1 I =■  1/K,  2 

IX  = I TMPt  I f 2)  - I TMPt  1 ) 

I Y = I TMPt  !♦  3) -I  TMPt  If  1) 

L 1 = I TMPt  I) 

L2= I TMPt If  2) 

DO  41  L=L I ; L2 
I A=L- I TMPt I ) 

4 1 PTFtL) = I TMPt If  1 ) fl FIXt  FLOATt I Y) * FL  OAT  t I A) /FLOATt  IX) ) 

LPT-TTEt  2000) 

RE  TO KM 
FND 


•*«♦*»*«♦•*»*»«*»»«********< 


i***ii***»»****»*i 


» * 
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C SUBROUTINE  TO  DRAW  PART  OF  AN  UNDERLINE 

SUBROUTINE  DRAWC I TMP,  K) 

COMMON  /A/WC A) /B/TCC 2) 

INTEGER  I TMP( l ) 

CALL  DELPI  C(  2) 

CALL  BGNPI C<  2) 

CALL  LINI ( ITMP( 1), ITMP( 2>,0> 

CALL  LINI H< O/O, A) 

N = K-2 

IF<K.EQ.3)  GOTO  2 
DO  I 1=3, N, 2 

1 CALL  L I N I C I TM  P C I ) , I TM  P C I l ) ) 

2 CALL  ENDPIC 

TC(  1 > = C FLOAT ( I TMPC  N) > - W(  1 ) > /C  WC  2)  - WC  1 > > 

TCC  2)  = ( FLOAT  ( ITMPCN+1))-WC3))/CWCA)-WC3)) 
CALL  WRITOLC 201, 0, 0, TC, 2) 

RETURN 

END 

C SUBROUTINE  TO  ri ST LAY  DATA 

SUBROUTINE  DISPLAC I A, L, H, KEY, NPI C) 

CIMENSI ON  I AC  I ) 

INTEGER  L, H 
CALL  BGNPI CCNFIC) 

CALL  LINI (L, I ACL), 0) 

DO  1 I = L,  H, KEY 
1 CALL  LINI C I,  1 AC  I ) ) 

CALL  FNDPI C 

RETURN 

END 

C**********************************************> 

C SUBROUTINE  TO  WINDOW  THE  DATA 


SUBROUTINE  WINCOWC I A, L, H ) 
COMMON  / A/W ( A ) 

INTEGER  IAC l),YL,YH,L,H 
DATA  I UAL/ 2 00/ 

YL= I AC  L ) 

YH=I ACL) 

CO  I K = L,  H 

IFC I ACK)  .LT.YL)  YL= I AC  K ) 

IFC  I ACK)  .GT.YH)  YH  = I AC  K) 

CONTINUE 

VIC  1 ) =L 

VIC  2)  =H 

WC  3>=YL-2 00. 

WC  A)  =YH 
CALL  WINDWCW) 

CALL  BGNPI CC 99) 

CALL  LINECWC  1 ) ,WC 3) , 0) 
CALL  LINI RC  0, I UAL) 

CALL  LINECWC I), 0,0) 

CALL  LINECWC 2), 0) 

CALL  ENDPIC 

RETURN 

END 
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WISp*<«.s»Sst(l0ALm^«"eal 

cum  TO  DD.Q  


c*  * * 

c 

c*** 


*****«*  + « + **J***t***l 

SUBROUTINE  TO  TABULATE  COEFF I Cl  ENTS  FOR  THE  FILTER 


rt'L-HOUTI  NE  HAMTAP 
COMMON  /C/TAPL ( 200) 

LATA  P/3.  I A 1 5926536/ 

A= - f 00. 

DO  1 K*  1*200 

TAL-LC  K > 0.  54*0. 46*  COSC  P*  A/200.  ) 

A=  A+  2 

RETURN 


ENL 


C * * * 
0*  * * 

c 

C *** 


****«  + «+  **  + + ***  + ++  +**  + + * 
*+****%.+  *«*+  + * + ***%  + + + + * + * 


STAHT  RAMP  FOR  FILTER 


• ♦**•♦*****♦*♦*** 


S'.’F  ROUT  I NE  HAMPGN  ( PSTH) 
INTEGER  BSTR(l) 

LO  1 1=1* 100 
AI  = 1 

P STHC I ) =BSTR(  10  1) 

RETURN 

END 


c * * * 

c*  * * 

c 

C*** 


c * * * 
c*  * * 

c 

c * * * 


***..♦..♦.*.**♦.*.**•****♦** 

Er.r  RAUF  FOR  FILTER 

*♦*»***♦*»*♦***.**♦*♦»**»*♦***.*.*. 

SUL P C U T I NE  HAMENC C E FTP) 

INTGER  ESTRC ! ) 

DO  1 1-1*1 00 

LSTRC I + 1 ) = C 1 00. -FLOATC  I ) ) *ESTHC  1 ) / 1 00. 

CONTINUE 

RETURN 

END 

****.***«**«***.*****»******* 
**.**«**.**.»«*'*»***:***'»** 

COPY  ROUTINE  FOR  FILTER 
****♦>.**»*♦.********»*.*♦* 

SUBFOUTI NE  H AM  CCP( TTS*  EE  TR) 

INTEGER  FT S C I ) * ESTRC  1 ) 

CO  1 1=1*1 00 
ESTRC  I 1 ) = PT  S ( I ) 

RETURN 
END 


£♦******.««***.*♦**.♦*♦*. *.****.i 

C*,********<,**  + *********  + 3 

C SUEPCUTINF  ”r0  FILTER  THE  UNDERLINE.  IN  PRACTICE  THIS  WAS 
C REPLACED  IY  AN  ECU I VAL ENT  AS SEMILEP  ROUTINE. 


C***************1******3 


SUBROUTINE  HAMFILC  ESTR) 
COMMON  / C/TAPL C 200) 

INTEGER  PSTHC 1 ) 

DO  1 1=1*2000 
SUM®  0 

DO  2 K=  1*200 

SUM  = SUM* LSTRC I +K- 1 ) *TAELCK> 
ESTRC I )= SUM/ 108 . 

DO  3 !®  1*2100 

BSTRC  220 1- I ) =BSTR<  2101-1  ) 

RETURN 
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Figure  10.  CORREK  subroutine 
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